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ARGON  PUFF  GAS  SOFT  X-RAY  LASER 

I.  Introduction 


The  successful  demonstration  of  soft  x-ray  lasers1-^  using  as  the 
lasing  medium  the  plasma  generated  from  high  power  laser-target  interaction 
has  provided  considerable  impetus  to  a  number  of  alternative  techniques  for 
producing  an  x-ray  laser.  One  such  scenario  involves  the  application  of 
pul3e  power  drivers  to  implode  gas  puff  plasmas.  Of  the  many  variations  in 
the  target  design  there  are  three  basic  concepts  currently  under 
investigation.  They  are:  (1)  a  single  cylindrical  annular  gas  puff,  (2)  a 
single  cylindrical  annular  gas  puff  imploding  on  a  central  core  plasma,  and 
(3)  similar  to  0)  and  (2)  except  that  the  imploded  plasmas  are  not  the 
la3ing  medium  but  a  source  of  intense  x-ray  emission  for  photopumping  a 

target  material.  Scenarios  (2)  and  (3)  are  discussed  by  the  authors 

and  by  T.  Hussey  elsewhere. Since  the  single  gas  puff  plasma 
is  central  to  all  three  scenarios,  it  is  essential  that  we 

understand  fully  its  radiative  properties  and  behavior.  Therefore,  the 
focus  of  this  investigation  will  be  confined  to  the  dynamics  of  a  single 
gas  puff  plasma.  The  intent  here  is  to  concentrate  on  the  radiation 

kinetics  while  characterising  the  implosion  dynamics  of  the  gas  puff  plasma 
by  means  of  the  simplest  plausible  description,  i.e.,  a  dynamic  pinch 

model.  Obviously,  such  a  model  ignores  much  of  the  observed  pathological 
behavior  of  these  plasmas,  but  it  does  provide  a  starting  point  for 


determining  whether  it  is  at  ail 
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density  conditions  favorable  to  x-ray  lasing 


The  influence  of  the 


development  and  growth  of  plasma  instabilities  and  their  impact  on  the 
uniformity  and  homogeneity  of  the  evolving  plasma  is  a  critical  issue  ar.d 
is  currently  under  investigation.  Preliminary  results  based  on  2-D 

a 

Magneto-Radiation-Hydrodynamic  simulations  and  experiments,  employing 
plasma  erosion  switches,  suggest  a  significant  reduction  in  the  role  of 
instabilities  and  provides  encouragement  for  further  pursuit  of  this 
design. 

The  treatment  presented  here  is  based  on  a  simple  model,  hereafter 
referred  to  as  SIMPL3DE,  with  emphasis  on  the  atomic  kinetics  and  radiation 
dynamics  self-consister.tly  coupled  to  a  dynamic  pinch  model  in  order  to 

obtain  information  on  the  plasma  environment  a3  it  evolves  under  the 
influence  of  the  current  discharge.  This  provides  us  with  a  self- 

consistent  picture  of  temperature,  density,  size,  and  level  population 
during  the  plasma’s  evolution.  The  influence  of  radiation  transport  on  the 
level  inversion  and  gain  calculations  is  taken  into  account  using  Sobolev's 
method.9,10  In  either  a  stationary  medium  or  one  moving  without  velocity 
gradients,  the  final  emission  of  optically  thick  line  photons  or  continuum 
from  the  medium  generally  occurs  from  unit  optical  depth  from  the 
boundary.  Photons  originating  deep  in  the  interior  of  the  medium  undergo 
successive  scatterings  on  their  flight  to  the  surface  before  escaping.  In 
this  instance  the  emitted  radiation  reflects  nonlocal  conditions.  Cn  the 
other  hand,  if  the  medium  ha3  a  velocity  gradient  then  radiation 

originating  deep  in  the  interior  of  the  plasma  can  escape  the  medium 

directly  because  of  the  Doppler  effect  and  this  radiation  reflects  the 
local  interior  conditions  rather  than  conditions  at  other  points  in  the 
medium  and,  ir.  particular,  the  boundary.  This  is  the  essence  of  Sobolev's 
approximation.  A  mere  thcrougn  discussion  of  the  effect  of  radiation 
transport  on  the  gain  of  the  lasing  transitions  can  be  found  in  reference 

2 


7.  Also,  since  the  experiments  on  the  GAMBLE  II  facility  at  NHL  are  being 
done  with  an  argon  gas  puff  plasma,  the  theory  and  analysis  presented  here 
is  for  argon. 

So  far  the  discussion  has  been  general  focusing  on  the  plasma  and  its 
properties.  However,  in  particular,  we  will  investigate  the  feasibility  of 
creating  a  population  inversion  in  the  3p  levels  of  r.eonlike  argor.  due  to 
electron  impact  collisions!  excitation  from  the  2p^  ground  state  and 
estimate  the  gain  coefficients  in  the  n-3,  An-0  lasing  transitions. 

II.  Physical  Model 

The  simplest  description  of  an  imploding  I-pir.ch  plasma  is  probably 
the  3er.r.ett  Pinch  equilibrium  model.  This  model  is  based  on  an  equilibrium 
balance  between  the  fluid  and  magnetic  pressures  in  combination  with  an 
equilibrium  balance  between  the  sources  and  sinks  of  energy.  The  classical 
Sennet  pinch  ignores  radiation  and  hence  ignores  excitaticn/ionization 
energy  (chemical  potential)  and  energy  lost  by  radiation  (radiation 
cooling).  An  obvious  extension  to  the  Ser.r.ett  pinch  equilibrium  model  is 
the  inclusion  of  the  flow  parameters  describing  the  temporal  evolution  of 
the  plasma.  That  i3,  maintain  the  simple  philosophy  of  the  Ser.r.ett  pinch 
but  allow  the  plasma  to  radiate  and  evolve  ir.  time.  In  essence,  this  i3 
our  philosophy  -  a  radiating  dynamic  Ser.r.ett  pinch;  3IMPLCDE.  Also,  like 
Shearer,'1  who  included  a  radiation  cooling  term  in  the  form  of 
3remsstrahiur.g  losses  from  a  pure  hydrogen  plasma,  we  include  radiation 
cooling  but  in  a  much  more  extensive  fashion.  The  radiation  cooling  term 
ir.  our  model  includes  contributions  from  free-free,  free-bour.d,  ar.d  bound- 
bound  transitions  ar.d  is  determined  from  a  r.or.-lTS  collisior.al-radi ative 

dynamics . 1 2 
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model  of  the  level 


The  SIMPLGDE  model  describes  the  radial  implosion  of  a  cylindrical 
annular  gas  puff  plasma  of  uniform  density  carrying  a  uniform  current  in 
the  axial  direction.  Only  radial  motion  is  considered,  i.e.,  there  is  no 
axial  structure  and  the  plasma  is  always  uniform  in  this  direction.  The 
radial  motion  is  determined  from  the  force  equation,  vis. 


r2 

d  r  i 

m  2  *  (P  ~  2  o 

dt  2irr  c 


where  r  is  the  radial  distance  measured  from  the  symmetry  axis,  P  is  the 

fluid  pressure,  A  is  the  area  over  which  the  force  is  exerted,  i.e., 

2  2  2  ? 

2 mrl,  and  I  /2irr  c  13  the  magnetic  pressure,  i.e.,  3  /St.  The  thermal 

energy,  S.»3/2  (1+Z)n.kT  +  Eu  ,  varies  in  time  as 

th  i  p 


P  V 
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where  Zu  is  the  sum  of  ionization  energies  and  is  loosely  -eferred  to  as 
? 

2  2 

chemical  potential.  A„urr  »  t(P^-S^)  where  represents  the  outer 

(inner)  radius  of  the  annular  plasma,  7»  2*rirr,  M»(m/m4  )/-rr“Z  where  n  is 
the  plasma  mass,  m^  is  the  atomic  mass,  n4  «  m/m4 ,  and  the  thermal  pressure 
is  N((+2)kT.  Finally,  1  i3  the  length  of  plasma,  3  is  the  charge  state,  z 


is  the  classical  resistivity,  and  ?ra(4  is  the  power  radiated  per  unit 
volume  V.  The  remaining  sysmools  have  their  usual  meaning.  The  first  term 
on  the  P.HS  of  Eq.  '2  represents  the  work  done  in  compressing  the  plasma, 
the  second  expression  is  the  joule  r.eating  source  term,  and  the  third  term 
is  the  power  radiated,  i.e.,  radiative  cooling. 

The  coilisionai-radiative  model  describing  tr.e  level  dynamics  contains 
an  extensive  number  of  levels  in  the  ;<-  and  ’--shells  with  particular 


—  to  _ 


h  cor.t 


represented  by  Voigt  functions  and  include  natural  and  Doppler 


broadening.  The  transport  of  radiation  employs  the  Sobolev  escape  model. 

The  " x-ray”  lasing  scheme  considered  here  takes  advantage  of  the  large 
monopole  excitation  rate  from  the  ground  state  to  the  excited  3p  state  of 
the  neonlike  ionization  stage.  Although  we  are  focusing  on  a  collisional 
excitation  scheme,  all  other  scenarios,  such  as  recombination  lasers,  can 
be  investigated  with  this  model.  However,  by  controlling  the  implosion  we 
hope  to  avoid  burning  through  the  neonlike  stage  and  instead  create 
conditions  for  producing  a  stable  abundance  of  neon! ike  argon.  A 
simplified  energy  level  diagram  for  neonlike  argon  is  shown  in  Fig.  1  .  A 
population  inversion  can  occur  in  the  3d  and/or  3p  levels  leading  to  gain 
and  lasing  in  the  3d~3p  and  3p~3s  transitions.  This  is  indicated 
specifically  for  the  3p~3s  transition  with  the  label  KU311A). 

III.  Results  and  Discussion 


Preliminary  estimates  for  the  neonlike  fractional  abundance  and  gain 
coefficient  ir.  a  stationary  environment  were  obtained  from  the  coilisienal- 
radiacive  model  for  prescribed  values  of  temperature,  density,  velocity, 
and  3iza.  These  estimates  crovide  a  measure  of  the  oarameter  soace  over 


which  gain  can  be  expect 

ed 

a3  well  as  provid 

ir.g  guidance  for 

simulations.  In  Fig.  2 

the 

neonlike  ground 

state  fractional 

abundance  is 

presented  as  a  function 

of 

temoerature  for 

several  ior.  der.s 

ities  in  the 

or 


absence  of  opacity  effects,  i.e.,  the  optically  thin  case, 
temperatures  from  about  30  to  7D  eV,  for  densities  typical  of  imploding  gas 
puff  plasmas,  a  significant  amount  of  r.eonlixe  argor.  prevails.  The  gain 
coefficient  for  the  3o-3s  transition  is  sncwr.  in  Fix.  3.  for  an  ootioallv 


a  *  10  f  A (A) (M/T)  (N_~  — —  N, )  cm  . 

oac.  2  g. 


where  f-_-  is  the  absorption  oscillator  strength  of  the  line,  A  is  the 
wavelength  in  angstroms,  M  is  the  atomic  mass  of  the  radiating  ion,  T  is 
the  temperature  of  the  plasma  in  eV,  N2  (M« )  is  the  upper  (lower)  level 
population,  and  g^  is  the  statistical  weight  of  level  i.  The  gain 
coefficient  for  a  Doppler  broadened  line  is  directly  proportional  to  the 
wavelength  of  the  lasing  line  and  the  difference  between  the  upper  and 

lower  level  population  densities,  and  inversely  proportional  to  the  square 
root  of  the  temperature.  For  an  optically  thin  plasma,  gain  coefficients 
greater  than  -unity  exist  for  all  three  densities  over  a  broad  temperture 

_  1 

range  with  a  peak  gain  coefficient  of  about  30  cm  '  at  50  eV  for  an  ion 

density  of  1 .5x10 '^cm'^.  For  a  fixed  density  and  increasing  temperature 

the  number  of  necnlike  ions  decreases  due  to  increased  ionization  causing 
bum through;  for  decreasing  temperature  the  plasma  becomes  too  cold  to 
support  the  existence  of  neor.lixe  ions,  and  the  colli3ior.al  excitation 

of  *he  '3o  r.eenlike  1  eve  Is  .iron  drastically.  Similar!/,  for  a  fixed 


temperature  ar.c  increasing 


abur.dar.c 


ions  decreases  due  to  t 


aus 1  ns 


burntr.rougn;  another  way  o 


density  increases  for  a  : 


temperat 


plasma  ■ 


GAM5L 


differences  are  greatest  for  the  higher  densities  due  to  the  combined 
effects  of  collisions  and  the  radiation  field. 

The  combined  effects  of  velocity  and  opacity  can  but  be  understood  by 
first  considering  them  separately  and  then  as  a  composite.  In  Fig.  6  the 
results  of  several  calculations  are  presented  for  varied  conditions  while 
maintaining  the  total  ion  density  fixed  at  5x10  'cm  J.  The  optically  thin 
result  is  included  for  reference  and  comparison  purposes.  The  influence  of 
opacity  on  the  production  of  neonlike  ions  is  3hown  on  the  v=0,  diameter  » 
0.9mm  curve.  In  comparison  with  the  optically  thin  case,  it  is  seen  that 
an  increase  in  plasma  size  manifests  itself  with  an  increase  in  opacity 
which  maintains  the  radiation  field  in  the  plasma  thereby  making  it  easier 
to  achieve  a  given  degree  of  ionization,  in  this  instance  the  neonlike 
stage,  at  a  lower  electron  temperature  than  the  optically  thin  case. 
Therefore,  for  a  fixed  electron  temperature  opacity  effects  will  enhance 


the  plasmas’  degree  of  ionization  above  that  achieved  in  a  purely 
coilisioral  optically  thin  plasma.  Hence,  the  rapid  decrease  in  neonlike 
ions  when  the  plasma  is  opaque.  With  increasing  temperature  coiiisional 
ionization  will  reduce  the  number  of  neonlike  ions  in  both  cases.  The 
influence  of  velocity  on  the  neonlike  ion  fraction  can  be  understood  in  the 
following  way.  In  plasmas  where  the  directed  motional  velocity  exceeds  the 
thermal  velocity  the  escape  probability  is  enhanced.  In  essence,  when  the 
line  Quanta  innings  on  a  region  of  cias.na  where  the  local  velocity  has 


associated 


ve.ocity 
•y  plasma 


e  ve.oc 

ro-oya 


he  enzire  plasma 
reduce  the  effective  opacity 
the  plasma  effectively  thin 
he  effects  of 


Finally,  to  explore  the  effects  of  mixing  elements  we  have  included  the 


results  of  calculations  where  lO1^  neon  ions  per  cm^  are  mixed  with  5x1 01^ 
argon  ions  per  cm’.  These  results  are  shown  also  in  Fig.  6  with  the 
consequence  of  further  reducing  the  number  of  neonlike  argon  ions  due  to 
increased  ccllisional  effects  induced  by  the  higher  density. 

The  gain  calculations  for  the  conditions  depicted  in  Fig.  6  are  shown 
In  Fig.  7.  Note  the  inclusion  of  an  additional  gain  curve  for  a  velocity 
of  1x10 'cm/sec  and  a  diameter  of  d-0.9mm.  For  the  conditions  described  in 

Fig.  7,  gain  coefficients  can  be  achieved  in  excess  of  unity  in  the  3p~3s 

lasing  transition  at  434A  for  peak  implosion  velocities  in  excess  of 
1x107cm/3ec.  It  is  clear  from  these  single  plasma  calculations  that  for 
prescribed  values  of  temperature,  density,  velocity  an  size,  that  are 
representative  of  gas  puff  implosions,  it  should  be  possible  to  create  a 
population  inversion  and  achieve  gain  coefficients  in  excess  of  'unity. 

We  will  now  investigate  whether  these  conclusions  will  prevail  in  a 

dynamic  environment  such  as  that  generated  by  the  SIMPL3DE  model.  For 

illustrative  purposes,  we  present  the  results  of  calculations  for  a 
reference  case  of  a  a  cm  long,  35  ugm/cm  argon  gas  puff  plasma  distributed 
uniformly  between  the  outer  and  inner  radius  of  1.55  and  0.95  cm, 
respectively.  The  driving  current  waveform  typical  of  the  GAMBLE  II 


generator,  with 

the  Plasma  Erosion  Opening 

owl  zee , 

i3 

3howr.  in  Fig.  S. 

Peak  current  of 

3.75x10^  amps  i3  attained  i 

n  scout 

50 

nsec,  decaying  to 

about  ox*  l3  amps 

in  another  70  r.3ec  and  then 

falling 

precipitously  to  zero 

in  ICO  nsec.  The  temporal  behavior  of  the  radiis  is  shown  in  Fig.  9.  The 
inner  radius  collapses  and  stagnates  on  axis  while  the  outer  radius 
continues  inward  until  the  back  pressure  is  sufficient  to  impede  the 
forward  metier,  and  bounces  outward.  The  final  pinch  radius  i3  about  1/10 
the  initial  racius  and  is  in  geed  agreement  with  the  oulk  of  experimental 
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data  accumulated  over  the  years  from  a  variety  of  generators  and  plasma 
loads.  The  variation  of  velocity  as  a  function  of  time,  shown  in  Fig.  10, 
reaches  a  peak  value  of  1.1x10^  cm/sec  in  about  160  nsec  which  is  well 
after  peak  current.  This  behavior  Is  characteristic  of  the  GAMBLE  II 
generator  and  is  especially  true  of  driving  currents  with  sharp 

risetimes.  The  implosion  phenomenology  suggests  that  the  plasma  heats  up 
and  percolates  for  a  time  and  then  eventually  coasts  inward.  The  temporal 
variations  of  temperature  and  density  are  shown  in  Figs.  11  and  12, 

respectively.  The  ion  density  peaks  at  the  pinch  and  reaches  a  value  of 

2x1 01®  ions/cm^  while  the  peak  temperature  occurs  some  15  nsec  earlier  and 
reaches  a  value  of  about  165  ev.  The  total  radiative  yield  which 
essentially  comes  from  the  L-shell,  is  roughly  1  kilojoule  and  exhibits  a 
pulse  duration  of  about  25  nsec  as  shown  in  Fig.  13.  This  result  is  in 
reasonably  good  agreement  with  the  experimental  observations  from  GAMBLE 


A  sample  of  the  gain  coefficient  for  the  43a  A  line  is  shown  as  a 
function  of  time  in  Fig.  14  for  the  illustrative  case  represented  in  Figs. 
3-12.  The  gain  coefficient  was  greater  than  unity  for  a  long  time  reaching 
a  peak  value  of  about  4  cm”1  late  in  the  implosion.  The  values  obtained 
for  the  gain  coefficient  are  probably  reasonable  around  peak  compression 


but  become  less  reliable  after  the  bounce  because 


:he  lack  of  an 


adequate  physics  description  of  this  late  stage.  However,  it  is 
encouraging  that  conditions  prevail  for  producing  measurable  gains  over  a 
4  cm  length  of  plasma  assuming,  of  course,  stability  of  the  column.  Work 
is  currently  in  progress  using  mere  sophisticated  mcdei3  to  assess  the 
validity  of  our  findings  here. 

Finally,  a  series  of  simulations  were  performed  to  determine  the 
lasma  parameters  at  maxi 
Dr  respond  in?  to  Hr-*'.  5: 


temperature,  at  maximum  gain,  are  shown  as  a  function  of  M/1  in  Figs.  *5 
and  15,  respectively.  The  gain  coefficient  as  a  function  of  M/1  is  shown 
in  Fig.  17.  For  Ar- 0.60,  the  gain  coefficient  peaks  at  M/l«60ugm/cm  and 
has  a  value  of  about  4  cm-1. 


IV.  Summary 


It  has  been  theoretically  demonstrated  that  it  is  possible  to  create 
conditions  favorable  to  population  inversion  and  gain  coefficients  in 
excess  of  -unity  for  a  variety  of  conditions  by  imploding  a  cylindrical 
annular  argon  gas  puff  plasma.  Gain  coefficients  of  4  cm"1  have  been 
calculated  for  the  3p-3s  lasing  transition  in  neonlike  argon  at  43U  a. 
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